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Abstract 

We study the propagator of a colored scalar particle in the background of a non-abelian 
gauge field using the worldline formalism. It is obtained by considering the open worldline of 
a scalar particle with extra degrees of freedom needed to take into account the color charge 
of the particle, which we choose to be in the fundamental representation of the gauge group. 
Specializing the external gauge field to be given by a sum of plane waves, i.e. a sum of external 
gluons, we produce a master formula for the scalar propagator with an arbitrary number of 
gluons directly attached to the scalar line, akin to similar formulas derived in the literature for 
the case of the scalar particle performing a loop. Our worldline description produces at the 
same time the situation in which the particle has a color charge given by an arbitrarily chosen 
symmetric or antisymmetric tensor product of the fundamental. 


1 Introduction 

The worldline representation of effective actions has seen a great deal of activity in the past 
twenty years, starting with the work of Strassler [T], who rederived Bern-Kosower master for¬ 
mulas lai directly from point particle path integrals, see ref. for a review. Since then many 
extensions and applications of the worldline formalism have been considered: multi-loop compu¬ 
tations |S], nonperturbative worldline methods EIITIIH], the numerical worldline approach to the 
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Casimir effect [5] , the worldline formalism in curved spacetime nnumna, photon-graviton mix¬ 
ing at one loop m, higher-spin field theory [HIIIIIII], the worldgraph approach to QCD 
as well as the heat kernel expansion |18l [191 EQ] , applications to noncommutative QFT [SUES, 
to the Standard Model physics [53] and its grand-unified extensions [53], just to name a few. 

Unlike effective actions, the worldline representation of dressed propagators is still a relatively 
unexplored land, though a worldline representation for the dressed propagator of a scalar field 
coupled to electromagnetism (scalar QED) was proposed long ago by Feynman j^. It consists 
of a worldline path integral where the coordinate paths have the topology of a line, and come 
with Dirichlet boundary conditions. The problem has then been reconsidered in ref. |26| along 
the lines of II1I1I3]. More recently, the dressed propagator in a scalar field theory has been 
studied with worldline methods to address the summation of ladder and crossed-ladder diagrams 
and analyze the emergence of bound states m- A full worldline description of dressed QFT 
propagators would be quite welcome in gauge theories as well. It may allow to address several 
different issues, providing a systematic way of computing scattering amplitudes that could be 
beneficial both at the perturbative and nonperturbative levels. At the perturbative level it may 
improve on the efficiency of perturbative calculations and give perhaps a better understanding 
of color/kinematics dualities |5S]. At the nonperturbative level it might be useful to address the 
emergence of bound states, or to study, for instance, the covariance of the Green’s functions under 
a change in the gauge parameter that gives rise to the the so-called Landau-Khalatnikov-Fradkin 
transformations |291 150] . Indeed the worldline formalism has been used recently to extend these 
transformations, originally derived for the scalar propagator, to an arbitrary n-point function in 
scalar QED m- The search of computational methods for tree-level amplitudes in gauge theories 
has been quite an active subject in the past decade. In particular, in scalar QCD, techniques 
based on recursive relations have been successfully found: see e.g. refs. |32l|33llSll35]. 

Here we consider the worldline approach to the propagator for a scalar field coupled to an ex¬ 
ternal non-abelian gauge field. Similarly to previous worldline treatments of non-abelian effective 
actions [331137], we obtain the path ordering—needed for the gauge covariance of the worldline 
path integral—through the quantization of suitable auxiliary variables [3311331130]. Their useful¬ 
ness resides in the fact that they allow one to get rid of the explicit path ordering prescription. 
Computationally it is a great advantage, analogous to the replacement of the path-ordered spin 
factors, present in Feynman’s original proposal for describing a spin- 1/2 particle |25| . with a 
Grassmann path integral over fermionic coordinates mi- However, the auxiliary variables thus 
introduced must be constrained in order to produce an irreducible representation (irrep) of the 
gauge group. One possibility is to couple them to a worldline U(l) gauge field, which allows for 
a Chern-Simons term whose coupling constant is chosen to project onto the desired irrep. This 
is the method already seen at work in the worldline description of p-forms El US], later ex¬ 
tended to the treatment of color charges. Furthermore it allows one to use commuting auxiliary 
variables which, in absence of the Chern-Simons term, would generate an infinite-dimensional 
color space. 

Below we show how to employ commuting auxiliary variables with a Chern-Simons term 
to study the propagator of a scalar particle in the background of a non-abelian gauge field by 
representing it as a suitable worldline path integral. The non-abelian charge (i.e. color) of the 
scalar particle can be arbitrary, though for simplicity we choose it to be in the fundamental rep¬ 
resentation. The propagator dressed by the external gauge field is (background) gauge covariant. 
One may then specialize the external gauge field to be given by a sum of plane waves (external 
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gluons) and perform the path integral. This produces a master formula generating the propa¬ 
gator with an arbitrary number of gluons directly attached to it. At fixed nnmber n of gluons, 
we call the latter “partial n-gluon scalar propagator”: indeed it consists of the scalar propagator 
with n gluons directly attached to it, while the glnon self-interactions are exclnded. As snch it 
is a gange dependent object. However this dressed propagator is valid off-shell, and it can be 
nsed as a bnilding block for higher-loop amplitudes. For example, being valid off-shell, one conld 
nse it to constrnct ladder diagrams with gluon rungs, and study the emergence of bonnd states, 
as done in ref. m for the pnrely scalar case. This would still be a gauge-dependent ladder, 
but in the bound-state stndies maintaining gauge invariance remains a difficult open problem. 
One may also try to nse it on-shell and stndy tree-level amplitudes, but then the addition of 
the redncible diagrams with three- and four-gluon vertices will be essential to achieve gauge 
invariance. We do not wish to discuss this final issue here, though one may hope that a kind of 
“tree replacement rnles”—such as those used in similar one-loop master formulas to generate the 
missing one-particle-redncible terms (see for example ref. |1])— might work here as well. Alter¬ 
natively, one may combine the present treatment of color charges with the worldgraph approach 
of ref. [T7], which is able to generate particle redncible graphs. We leave this to fnture analysis. 

The main point for ns in the present paper is to nse and exemplify a novel representation of 
the color degree of freedom. The partial n-gluon scalar propagator is the most natural starting 
point for testing its nsefulness in QCD applications. 

A final comment on the Chern-Simons conpling: fine tuning such coupling allows one to 
give the scalar particle a color charge corresponding to any arbitrarily chosen symmetric tensor 
product of the fundamental representation, rather than to the fundamental itself. Replacing 
the commuting variables with anticommuting ones would also give a similar construction, bnt 
with the non-abelian charge sitting in an arbitrarily chosen antisymmetric tensor prodnct of 
the fundamental representation, including the fnndamental itself. The choice is again made by 
selecting an appropriate Chern-Simons coupling. 


2 Tree-level amplitudes in scalar QED 


As a warm-up exercise, and in order to fix our notation, we review the compntation of tree-level 
amplitndes in scalar QED from a point particle (worldline) path integral on the line, i.e. with 
Dirichlet bonndary conditions. 

It is well known, since the seminal work of Feynman 125], that the propagator for a massive 
charged scalar field coupled to electromagnetism can be obtained from a worldline path integral 





DxDe 
Vol Gauge 


g-Sp,e;A] ^ 


( 2 . 1 ) 


where the particle action in Euclidean signatnre, is given by (space-time indices are left implicit 
where not reqnired) 


S[x,e;A\ = J dr (^^x^ + - iqx ■ A{x)'^ , (2.2) 

with e being the einbein, i.e. the gauge field for the one-dimensional diffeomorphisms. After the 
gauge fixing e = 2T, the path integral over e reduces to a numerical integral over inequivalent 
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constant gauge configurations labeled by the proper time T, i.e. 


pOO px{l)—x / \ 

{(l){x)^{x'))^= dTe-^^" 25a. e-/o. (2.3) 

Jo Jx{0)—x' 

In particular, treating the external electromagnetic potential as a perturbation, the expres¬ 
sion (2.31 can be shown to contain the sum of an infinite number of tree-level Feynman diagrams 
with an incoming scalar particle in x', an outgoing scalar particle in x, and an arbitrary number 
of photons. Specifically, in order to extract the amplitude with n photons, one first writes the 
potential as a sum of the n photons with polarizations and momenta ki, 


A,{x{T)) = Y^e,,i , 


(2.4) 


1=1 


and expands the exponential involving the potential, then extracts the amplitude as the term 
in (2.3) that is multilinear in all the different polarizations e^s. The amplitude thus reads 

pOO ^ 

A{x',x]ei,ki ,. ..SnAn) = {iqT / JT TT / dri 

Jo Jo 


rx{l)—x 


/ x{0)—x' 




m.l. 


(2.5) 


where the x ■ e terms have been re-exponentiated and xi := x { ti ), whereas “m.l.” stands for 
multilinear. We may now perform the path integral by splitting the generic path x { t ) into 
the background, Xhgir) = x' + {x — x ') t , satisfying the boundary conditions, and quantum 
fluctuations 2 /('r) with vanishing boundary conditions. We thus get 


r 

A{x', X ] £1, fci, . . . £„, kn ) = (ig)” / 

Jo 


dT 


{AttT)- 




n / 


1^1' 


X e 


J2i [iki-{x'+Ti{x-x'))+ei-(x-x')\ / J2i{iki-yi+eryi) 


( 2 . 6 ) 


where the expectation value ( • • •) is taken with respect to the free Gaussian path integral 

Iy{o)=o e“3T /o . Expression (2.6) is thus the expectation value of the product of photon 
vertex operators 


VaIeA] = ebk-ir{x-x')+y)+e-y] 

Jo 


From the free path integral one may obtain the worldline propagator 

{y^{r)y^^'{T')) = -2TS^^'AiT,T') 
A(r,r') = tt ' + - ^(r-tr'), 


(2.7) 

( 2 . 8 ) 
(2.9) 
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and eq. (2.6) can thus be written as 


A{x',x;ei,ki,...en,kn) = 

n 

X exp I ^ \iki ■ {x' + Ti {x 
1=1 




dT 




lo (iirT) 2 
- x')) + £i- {X- x')] 



+ (fci • ki'Aii> 


i2ei ■ ki>*Aiii 


£i 



( 2 . 10 ) 


1,1’ 

where A/;/ := whereas left and right bullets indicate derivatives with respect to t; and 

Ti' respectively. Therefore, 


'All, = Ti' - din, - n) 

•A'li, = 1 - Sin -n,), 


( 2 . 11 ) 

( 2 . 12 ) 


where 9{x) is the Heaviside function. Notice that—although products of delta functions appear 
in the full expansion of the exponential in the expression (2.10)—the multilinear part of the 


expansion, that yields the amplitude, does not involve singular, or ill-defined, terms. 

In order to get the amplitude fully in momentum space one may Fourier transform the 


expression (2.10), and get 


Aip',p;£i,ki,...£n,kn) = J j h, ■ ■ ■ £n, kn) 

= J d^X+ j d^x- A^P+P'A+ + Up-P'y--Aix+ - ^,X+ + ^-,£1^1,. . .SnAn) , ( 2 - 13 ) 

with X- := X — x' and x+ := The integral over the “center of mass” x+ yields the energy- 

momentum conservation delta-function { 211 )^SA) '^ki^ , whereas the integral over the 

“distance” x_ is Gaussian. Hence, after some simple manipulations, the amplitude reduces to 

_ pOO ^ pi 

A{p\p]ei,ki,...€n,kn) = (iq)'^ / dT TT / dn 

Jo 

n 

X exp |r(p - p') • ^{-km + isi) TT^[ki- ki>Ai_i> - 2isi ■ A://A/_// + Si • | 

i.v 


m.l. 

(2.14) 


with := i |t;—T;/ 1, i.e. the translation-invariant linear part of the worldline propagator Q, 
a fact that was already noted in ref. [26]. For simplicity, and for later convenience, in eq. (2.14) 


the overall energy-momentum conservation delta function has been stripped off, and the resulting 
amplitude has been indicated with A. 


Equation (2.141 is the Bern-Kosower-type formula for the tree-level amplitude with n photons 


and two scalars, originally worked out in ref. [2^ (see also ref. [3T]). Note that this amplitude is 
amputated on the external photon lines, and may be called the n-photon scalar propagator. It 
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gives the sum of all the QFT Feynman diagrams depicted in Fig.[2 which involve linear vertices 
and seagull vertices, and where the external photon momenta appear in all possible orderings. 
This is an advantage of the worldline formalism, that it combines all orderings into a single 
expression and that, for scalar QFT, generates the seagull vertices within the linear worldline 
vertex operator. It is easy to check that, for n = 2, it reproduces the sum of the three Feynman 



Figure 1: Diagram structure for the n-photon amplitude 


diagrams responsible for the scalar QED Compton scattering. Notice that the integrand in 
eq. (2.141 is made of a worldline translation-invariant double sum, and a simple sum that is not 
translation-invariant. However it becomes translational invariant on the shell of the scalar lines: 
indeed, imposing the translation Sti = A on such sum, and using momentum conservation, it 
gives A(p — p') ■ ~ P^)) which vanishes on-shell. 


3 Dressed scalar propagator in a non-abelian background 

A straightforward way to generalize the previous results to the case of the coupling of a scalar 
field to a non-abelian gauge field = IF“(x)T“, with T“ belonging to the Lie algebra g 

of a gauge group G, is to use a path-ordering prescription to guarantee the gauge covariance 
of the path integral. Let us thus take a colored scalar field with an index in an irreducible 
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representation of a gauge group, which for definiteness we assume to be G = SU{N). The scalar 
field propagator in the presence of the non-abelian gauge field can thus be written as 




c(l)=a: 


DxDe 


,:(0)=x' Vol Gauge 


(^g^-S[x,e-,W] 



(3.1) 


where the particle action in Euclidean signature is given by 

S[x,e\W]= j dri^^x^+ - igx ■W{x)^ , (3.2) 

and with V denoting the path ordering. Although this expression contains implicitly the multi¬ 
gluon scalar propagator, it may be convenient to investigate an alternative approach that im¬ 
plements the path ordering using some auxiliary variables. For the case of one-loop gluon am¬ 
plitudes, generated by a scalar field or by a Dirac field, this procedure was proposed in ref. m 
and more recently reconsidered in refs. |36ll37]. The main advantages of such a procedure are: 
(i) automatic implementation of path ordering, through the quantization of the auxiliary fields, 
that allows one to simplify the perturbative calculation, (ii) Simple way to extend our formulas 
to the case of a scalar field in a generic (anti)symmetric tensor product representation of the 
gauge group. 

For definiteness, let be in the fundamental representation, i.e. (T“)q,“ , a = 
then the totally (anti)symmetric tensor (j)ai...ar transforms in an irreducible (anti)symmetric 
product of r fundamental representations. Let us thus consider the addition of complex auxiliary 
variables, Ca and c“, that sit in the fundamental and antifundamental representations of the 
gauge group, respectively, and that are coupled to the gauge field potential 

/■l r 


/ dr c^c^-igx^W^c‘^{T‘^)Jcp 


- c“cq(1) . 


(3.3) 


JO 

Note the boundary term that allows us to set initial conditions on c and final conditions on c. 
The canonical quantization of such fields gives rise to creation and annihilation operators c and 
that generate a suitable Fock space. Using a coherent state basis like the one summarized in 
Appendix]^ it yields a scalar wave function in the coordinate u, the left eigenvalue on coherent 
states of the creation operator c^, that in turn can be expanded as a u-graded sum of fields 


(l){x,u) = (t){x) -f U°‘(l)c,{x) + 


(3.4) 


so that if the auxiliary variables are (anti)commuting, the tensors 4>ai...ar. sit in an (anti)- 
symmetrized tensor product of r fundamental representations. In the following we will concen¬ 
trate on the case of commuting auxiliary variables, as the anticommuting counterpart is just a 
straightforward generalization. Another necessary ingredient is thus a projector that allows to 


single out from eq. (3.41 only an irreducible tensor with r indices. This may be achieved by 


adding a C/(l) worldline gauge field to the action, along with a Chern-Simons term. The full 
locally symmetric action reads 


Si[. 


X, c, c: e, a 




-k 


m^e 


c“ 


(T“ 


/3 


C/3 


'(5r -k ia)Ca — isa^ — c“Ca(l) , 
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(3.5) 






where a(r) is the C/(l) gauge field and s is the Chern-Simons coupling, to be fixed shortly. Due 
to the presence of the auxiliary fields, the potential is no longer matrix-valued and no explicit 
path-ordering is thus needed. The propagator reads 




.x(i)=x,c(i)=u DxDWcDeDa 
Ix{o)—x'^c{o)—u' Vol Gauge 


^ — Si [x,c,c,e,a-,W] 


(3.6) 


where u'^ and Ua represent the initial and final color states of the scalar field. Thus, if the scalar 
field sits in the rank-r totally symmetric representation, the above field operators are (apart 
from numerical factors) given by 


(p{x, u) ~ (pai---ar{x)u°‘^ ' ’ ' 


^“i-“’-(x')n' •••< 


(3.7) 


Upon gauge-fixing the einbein e can be set to a modulus 2T as above, whereas the abelian 
gauge field can be set to a constant angle modulus 9 G [0, 27 r). The gauge-fixed path integral 
thus reads 


((/>(x,u)(/>(x',u'))^ 


^oo p2Tr in /•x{l)=x i.c(l)=u 

/ dT ^ / Dx / DcDc ^ 

Jo Jo Jx{0)—x^ J c(0)—u' 


(i)=fi 

(0)=t, 


(3.8) 


with 

S[x,c,c-2T,9,W]= dT(^^x^-tgxf^W;[c^{T<^)Jcp + c‘^{dr+t9)c^'^-c‘^Caiil) , (3.9) 

where the constant part —Tw? + isO was stripped off. Classically the equation of motion for a 
corresponds to the constraint c“Cq, — s = 0. In the canonical quantization approach this turns 
into an operatorial constraint to be imposed d la Dirac-Gupta-Bleuler on the wave function. As 
shown in the Appendix [^we choose to identify the path integral expression c^“Cq, as the one 
representing the normal-ordered hamiltonian cI“cq. Therefore the quantum constraint reads 

(/)(a;,c) = 0 , (3.10) 


and thus, in order to single out the tensor with r indices, we identify s = r, and set the path 
integral normalization to (see eq. ( A.24[)) 


pc{l)—U ( K 

/ £)c£)c e’-^o 

J c(0)— 


.( 1 ) 


(3.11) 


with u ■ u' = In the following we find it convenient to twist the auxiliary fields in order 

to absorb the 9 term, i.e. c(t) —>■ c(T)e“®®'^ and c(t) —>■ c(r)e*®'^, so that the previous expression 
gets replaced by 


^c(l)=e 

c{0)—u' 


DcDc e -^0 


— dr c“Cq,+c“Cq, (1) _ ^^u-u' 


(3.12) 
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and (3.81 reduces to 


{(t){x,u)(l){x',u'))^ = dT e 


-Tm^ 


r2TT 


de 

27r 


Ar0 


l'x{l)—x pc{l)—e 

Dx 


/ x{0)—x' 


DcDc 


c{0)—u' 


X e 


- dT(^^x^-igx^W^ c“{T“)a^c^+c“cQ)+c"ca(l) 


Splitting as usual the fields into backgrounds and fiuctuations 

Ca(T) = + Kq,(t) , K„(0) = 0 

c“(t) =e-*®it“ + K“(T), k“(1) = 0, 

the kinetic action for the fluctuations fg drR'^ka can be inverted to give the propagator 

{Ka{T)Rf^{a)) = 5^e{T-a) . 


(3.13) 

(3.14) 

(3.15) 

(3.16) 


Eq. (3.131 gives the final form of our worldline representation for the scalar propagator 


dressed by a non-abelian gauge held. It is gauge-covariant, and we discuss the related Ward 
identities at the end of next section. For an arbitrary external gauge held configuration one is 
not able to perform the path integral exactly. On the other hand, one may use it for suitable 
perturbative expansion, or to get exact results for specific external held configurations. 

As our main application we specialize the external gauge potential to be given by a sum of 
plane waves. This allows for an exact path integral calculation, and produces a master formula 
for the scalar propagator with an arbitrary number n of external (amputated) gluons lines 
directly attached to it, generating what we have called the “partial n-gluon scalar propagator”. 
As already discussed, the final result will not be gauge invariant, but it may be used as a 
starting building block to get eventually gauge invariant amplitudes. Also note that our master 
formula is similar to the one-loop Bern-Kosower-type of master formulas used in scalar QCD, 
that contain all one-particle irreducible diagrams with a scalar loop and an arbitrary number 
of external (amputated) gluon lines, the only difference being that instead of a loop we have an 
open line for the scalar particle. 

Thus, in order to extract our partial n-gluon scalar propagator from the gauge-fixed worldline 


path integral (3.131 we write the potential as a sum of n gluons with polarizations ei, momenta 


ki and colors a;, i.e. 


W^{x) = 


(3.17) 


Z=1 


and this allows us to read off the gluon vertex operator from the interaction term, namely 


Vw[e,k,a] ■- ^[^k-ir{x-x')+g)+e-g] -c, ^ ^ 

Jo 


lin. e 

(3.18) 


which differs from the photon vertex operator of eq. (2.71 in the presence of the color factor, 
given by the color generator contracted with the auxiliary fields. The n-gluon two-scalar 
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term reads 


Aw{p,u;p',u';ei,ki,ai,...,£n,kn,an) = {ig)"^ j d^xd^x' ^ 

n „i 

n / dr. 


dT 


^-Tm'^-^{x-x' f I ^ Jre+. 


Jo (AnT) 2 Jo 

X ^ [kji , £^ , ^ , 


/=! ‘ 


(3.19) 


where the average is computed with the propagators (2.91 and (3.16). Due to the fact that the 
color part of the gluon vertex operator factors out, the integrand of the non-color part of the 
expression (3.19) is identical to its photon counterpart given in ref. (2.14). We thus obtain the 
following final result for the tree-level partial n-gluon scalar propagator 


{ P : P 7 ^ 5^17 7 7 • • ■ 7 7 7 ^n) — (^i?) 


dT / dn 


1=1' 


X exp |r(p - p') ■ -h i£i) +T'^[ki- ki,ti_i> - 2i£i ■ ki>ki_i> + ei ■ £i'ki-i'] | 

i Id' 


de 

27r 


irO-\-e 


+ S“Hri))(r“0c7/^ K + ^pAri)) ' '' 

+ K“"(T„))(T“")„/"(n^„ +K;3 „(t„))) , (3.20) 


(e-*®r 


whose corresponding QFT Feynman diagrams have the same structure as those depicted in 
Fig.[T] for the n-photon amplitude. Let us underline that, since no three- and four-gluon QCD 
vertices are involved in our worldline computation, expression (3.20) only represents the gluon- 
irreducible part of the full n-gluon two-scalar amplitude, i.e. the part of the amplitude made of 
Feynman diagrams that cannot be parted by cutting a gluon internal line. 

Before studying some special case let us consider the 9 integrals present in the previous 
Bern-Kosower-type formula. Let us define 


r2-K 


Fik,u-u'):= I ^ e^ke+e '‘>u-u' 


(3.21) 


with k ^ 1.. We can solve it by transforming it into a clockwise contour integral over the unit 
circle: let us define z := Hence, 


F{k, u- u') = 


dz 


1 


—27rfz 


Uu-u'f 

0, 


fc > 0 
fc < 0 


(3.22) 


Notice that, although the function e“'“ ^ presents an essential singularity at z —>■ oo, it is a 
holomorphic function in any bounded domain, and it can be written as a power series. Thus, 
the contour integral over the unit circle picks out the order k coefficient of such power series. 
Let us now single out a few examples in order to further clarify the results obtained. 
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3.1 Free scalar propagator 

The free scalar propagator, for a field in a totally symmetric rank-r representation, can obviously 
be obtained from the above formalism by considering zero external gluons (n = 0). This will 
help us to fix an overall prefactor. For n = 0, eq. (3.201 reduces to 

1 


AwiP, u; -p,u') = 


P‘‘ + 771 ' 


F{r,u ■ u') = 




r-i 


+ 777'^ 


(3.23) 


where 5aX- '^ is the identity in the totally symmetric rank-r representation. 

Equation (3.23) corresponds to the free propagator for the polarized scalar field, and one 
may obtain the free propagator for the unpolarized scalar field by stripping off the prefactor 

~—from (3.23). In fact note that 


U ^ ■■■u ’ Uo 






■Ur. 


l[.'t=Pr{u\u') = Pre^-^ , 


(3.24) 


where Pj. is the projector on the above irrep—basically it picks up the order r term in the Taylor 
expansion of the exponent. Furthermore, using the identity 


1 

r! 


N 




Q!=l 


2777 


'Oil 


■ Ua = 


1 


dr 


dr 


r! • • • dv°‘r dvp^ ■ ■ ■ dvp^ 


v—v—0 


'-'Cti-'-OCr ’ 

one can easily prove the following composition rule for the free propagator 

/. /. N 




Pj. 


a—1 

P.r 


2717 


p^ + 


\x', u'){x', u’' 


Pr 


p^ -\- 777 ' 


|x",77") 




\x\u") . 


(3.25) 


(3.26) 


In eq. (3.25) we used the sources Va and 7 ;“ to insert the terms Ua^ ■ ■ ■ • ■ ■ u^'' into the 

Gaussian integral. 

Of course the prefactor -- — is universal, i.e. it only depends on the two external 

scalar lines and not on the gluon vertices. Therefore the recipe to get the generic 77 -gluon two- 
scalar term from eq. (3.20), is to strip off such color prefactors and truncate (i.e. multiply by 
the inverse free propagators) the external scalar lines. 
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3.2 Partial n-gluon scalar propagator 

For the simplest case of one gluon, expression ( |3.20| reproduces the Feynman vertex linear in 
the gluon, whose color part is given by the expression 


where 


F^[u,u') = ^ 1 (F“z + S“(r)) + npir)) 

^ • u'y-^uT'^u' 


—2'Kiz z'' 1 


(r- 1)! 


— (T'^ ... 7 / 




(3.27) 


(3.28) 


is a generator of SU{N) for the rank-r totally symmetric representation. Above we used the 
fact that {kp{t)R°‘{t)) = ^6^ and that r“ is traceless. Thus, it reads 


pOO pi 

Aw{p,u;p',u';e,k,a) = ig dT + p ) + 2 p fer) 

JO JO 


X ie ■ {p — p') 




r! 


'('7^0 N 01---0r„/ ... 1 /' 

VJ (r)/ai...ar “/3i 


= *5 


is - {p- p') 


(p2 + m?){p''^ + m^) 


('T’O \Pl---dT 
(r) /ai.. .a; 


. . . u°‘^u'p^ ■ ■ ■ u'p^ 


(3.29) 


that—after truncating the external scalar lines and stripping off the factor -correctly 

reproduces the one-gluon two-scalar vertex 

^ 2 s,ig{p,a\p\y,s,k,a) = igis ■ {p - p'){T^r))ti::^y ■ (3-30) 

For the case of two external gluons the auxiliary field correlator involves two generators, and 
the color part reads 


F“i“=(u,m') = ^ 


dz 1 


(u“i0-k«;“i(ri)) (r“i)a/'(M/3i +KfiyTi)) 


—2Triz z'' 

(u“=Z-kK“"(T2)) (T“=)a/"(u/32 +KpyT2)) 


dz 

— 2TTiz 


^ uT“iu' uT“=u' -k 6 I(ti - T2) + 

+ ^ uT^^T^^u' 0{t2 - Ti) + - tr(r“ir“=)0(Ti - r2)0(r2 - n) ) . 


(3.31) 
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The last term obviously vanishes for all r’s, whereas the first-one is non-zero when r > 2. 
Furthermore, by introducing 1 = 0(ti — T 2 ) + 9{t2 — ti) in front of the first term, the latter 
reduces to 


1 


(r- 1)! 
1 


Sr>2{r - l){uu'y-^uT^^u' uT'^^u' + {uu'f-^uT'^^T'^^u' 6I(ti - T2 ) 


SrMr - l){uuy-^uT'^^u' uT^^u' + {uu'y-^uT^^-^T^^^u' 6 »(t2 - n) , (3.32) 


(r- 1)! L 

with an obvious notation for 5r>2- The expressions in the square brackets are just the products 
of the generators in the rank-r totally symmetric representation, decomposed in terms of the 
fundamental representation, i.e. 






+ 9 (r,-r.)(T«r“;,) 


/31.../3. 

Cfl...OLr 


(3.33) 


that, once again stripping off the prefactor , gives the correct color factors for the two-gluon 


term of the dressed scalar propagator. Finally, one can insert the latter into the formula (3.201 


and get the final result for the partial two-gluon scalar propagator: the integrals over ti, T 2 and T 
reproduce the correct form factors. In particular the term involving the expression £i -£2 5{ti—T2) 
corresponds to the Feynman diagram with the two-gluon “seagull” vertex. 

The partial n-gluon scalar propagator term turns out to be a straightforward generalization 
of the previous expressions. In particular, for an even number of gluons, the color factor displays 
the same features as the two-gluon one: One finds trace terms that vanish as they multiply a 
full product of theta-functions. All the other terms combine into time-ordered products of color 
generators in the rank-r totally symmetric representation, i.e. 




7;“i • . . 


73i 


X E ^(^-(1) - • • ■ ^(^-0-1) - r.in))iTpy^ ■ ■ ■ , (3.34) 


<y^Sn 


where Sn is the group of permutations of n elements. The same happens for an odd number of 
gluons. 


3.3 Gauge-covariance of the dressed propagator and Ward identities 

The worldline representation for the dressed propagator in scalar QCD, eq. ( |3.13[ ), comes in 
handy to show the gauge covariance of the propagator under finite transformations. 

Upon a finite gauge transformation U{x) = with e{x) = e“(a:)T“ and in the 

fundamental representation, the gauge field transforms as 


Wf.{x) 


U{x) 




C/t(x) , 


(3.35) 
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and the worldline lagrangian in (3.13) is gauge-invariant provided the auxiliary fields transform 


as 


c(t) —>■ U{x{t))c{t) 

c{t) —)■ c{t)U\x{t)) . (3.36) 

In particular, using that u' and u are the boundary values for c(r) and c(t), at r = 0 and r = 1 
respectively, the previous rule implies that 


< ^ iU{x'))a°‘u'^=:v'^ 

^ u^(U^{x))/=:v^ , (3.37) 


and 


4){x,v)4>{x'= {4>{x,u)(l){x',u')^ ^ 


(3.38) 


for the polarized scalar propagator. Using the transformation rules (3.37) and stripping off the 
color factors, we obtain 


((^5i...5,(a;)(^^i'"^'-(a;'))~ 

\ / w 

= {U{x)U,^^ ■ ■ ■ {UHx'))p/^ ■ ■ ■ iU\x^))p/’^ , (3.39) 

\ / w 

that is the correct transformation rule for the dressed propagator of a scalar field in the rank-r 
totally symmetric representation of SU{N). At the infinitesimal level, the previous expression 
yields a generating functional of Ward-identities for the above amplitudes, namely 

• (3.40) 

For example, at the leading order in the coupling constant g the latter gives an identity between 
the one-gluon two-scalar term and the free scalar propagator. For definiteness let us consider 
r = 1. In Fourier space, and using the truncated amplitudes, the previous functional identity 
reduces to the perturbative Ward identity 


0 = 5^^'> {k + p + p') A 2 s,ig{p,a;p', P; -ik, k,a) + ig{T°-)a^ {p'^ -p'^) 


(3.41) 


that does indeed vanish. Notice that in (3.41) the replacement e —>■ —ik implements the covariant 
derivative of the first term of (3.40). 
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4 Conclusions and outlook 


We presented a worldline model where a scalar relativistic particle is coupled to a non-abelian 
gauge field. The quantization of the model yields the scalar propagator dressed by the external 
gauge field. The path ordering—needed for the gauge-covariance of the model—is realized using 
(commuting) auxiliary fields that are coupled to the gauge potential. The addition of a worldline 
abelian gauge field along with a Chern-Simons term allows one to describe scalar particles that 
sit in an arbitrary rank-r symmetric representation of the gauge group. We have specialized the 
non-abelian gauge field to be given by a sum of plane waves, and obtained explicit expressions 
for the scalar propagator dressed by an arbitrary number n of external gluons directly attached 
to the scalar line, the partial n-gluon scalar propagator. This is, in general, a gauge-dependent 
object, as three- and four-gluon vertices are not included, but it is valid off-shell and it could be 
used as a building block for constructing amplitudes. The resulting expressions are factorized 
into a color part and a kinematic part: their integral representation, provided by the present 
worldline model, may nicely combine with integration-by-parts techniques, that were shown to 
simplify the calculation of QCD form factors gsmiiig. 

One future line of research for extending our work would be to find a simple way of adding 
the gluon self-interactions and obtain the complete tree-level amplitudes with 2 scalars and n 
gluons. One option would be to investigate a kind of “tree replacement rules” of the type used 
successfully in similar Bern-Kosower formulas, where the scalar particle performs a loop instead 
of a line. 

Also, the present approach can be generalized straightforwardly to scalar QCD, with scalar 
particles in a rank-r totally antisymmetric representation, by using anticommuting auxiliary 
fields. Another quite welcome extension would be the treatment of the quark propagator in 
QCD, that may be obtained by coupling the external non-abelian field to a spinning particle, 
i.e. by using a particle with worldline supersymmetry, where the quantization of the spinorial 
coordinates describe the spin degrees of freedom. An alternative approach, that needs path 
ordering, is the inclusion of a matrix-valued (spin-factor) potential—see e.g. refs. |25[ITT] . 

A possible application of the n-gluon master formula obtained above, would be the recon¬ 
struction, in any covariant gauge, of the scalar-gluon vertex which is known as the Ball-Chiu 
vertex |49]. In previous studies of form factors for QED and QCD—see refs. laillallll] —it was 
shown that the worldline formalism simplifies such calculations, and the results can be expressed 
in a much more compact way. Finally it would be nice to extend the present treatment by adding 
the coupling to external gravity, along the lines of what has been done for the one-loop effective 
actions in curved space nni nulla. 
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A Coherent state path integral normalization 


Let us consider a Hilbert space spanned by a complete set of harmonic oscillator states {|n)} 
with n = 0 ,l,-- - , 00 , and let and c be a pair of creation and annihilation operators acting on 
that Hilbert space and satisfying the canonical commutation relation [c, c^] = 1. These operators 
generate the harmonic oscillators states {|n)} in the usual way. Let us now define the ket and 
bra coherent states as right and left eigenstates of c and c^, respectively 

c\u) = u\u) , {u\c' = (ttlu , (A.l) 

where u and u are complex numbers. They can be constructed out of the Fock vacuum as 

|n)=e->), {u\ = {0\e^\ (A.2) 

and they form an overcomplete basis of the Hilbert space. Indeed their overlap reads 

{u\u') = . (A.3) 

After this introduction to bosonic coherent states we prove the following identity: 


{u\e- 


y) = 


(A.4) 


A simple way of proving it is to consider the wave function (ulip), corresponding to a generic 
state |'0) of the Hilbert space. On this wave function the creation and annihilation operators act 


as c 


f 


u and c ^ so that the operator e 


—iOc^c 


e aa is just seen to generate a finite 


scaling of the u coordinate, namely u —> e Then (A.41 follows immediately. An alternative 

(A.5) 


proof goes as follows. Let us start from the “normal order” rule 

m 

yer = Y,s{m,k){cy{c)'^, 

k^O 

where S{m, k) are the so-called “Stirling numbers of second kind” that are defined by 




i=0 


{k-iy 


(A.6) 


with iS'(n, 0) = Snfi and 5'(n, 1) = 1. Now from the definition and properties of coherent state 
and Eq. ( A.5[) one can write the left hand side of Eq. (A.4) as 




') = («l E 


n—Q 


= e““' ^ 


yey 


S'(n, k){uu'y 


n=0 




= 6“' f; y^zuuu') 

n' 


(A.7) 


n=0 


^These numbers appears often in reordering problems, as for example in the worldline two-loop computation of 
the Euler-Heisenberg effective lagrangian of scalar QED and spinor QED |48) . 
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where Tn{x) are the so-called Touchard polynomials. Such polynomials can be also obtained 
through the exponential generating function 






n—O 


(A.8) 


Therefore, the series in (A.7l yields the right hand side of Eq. (A.4|, which is thus proved. 


Expression (A.4| can thus be used to fix the normalization for the harmonic oscillator coher¬ 


ent state path integral. A free coherent state path integral can be built up by inserting spectral 
decompositions of the identity 


1 = 


dc dc 
2TTi 


e—|c)(c| 


into the expression (|A.3| for the scalar product. One gets 

1 — 1 ^ ^ C-kipk ^k—1^ 


{u\u') = e““ = 


n dckdck 

exp c„c„_ 


n —1 


k=l 


2Tri 


fe=i 


(A.9) 


(A.IO) 


where we defined cq := u' and c„ := u. Furthermore we define a time parameter 0 < r < 1 and 
identify c(rfe) := Ck and c{Tk) := Ck, with — Tk-i = =: e. Thus, in the large n limit we 

may identify the latter as the free coherent state path integral 


with 


{u\u') = e®“' 



DcDc , 


Sf[c,c] = f dr c(r)c(r) — cc(l) . 

Jo 


(A.ll) 


(A.12) 


We may thus split the paths into backgrounds, satisfying the free equations of motion c = c = 0 
with corresponding boundary conditions, and quantum fluctuations 


c(r) = u' + k{t) , k(0) = 0 

c(r) = M -I- R{t) , k(1) = 0 . 


(A.13) 

(A.14) 


Hence, 


^ k (1)=0 

{u\u') = e““' = e®“' / DRDk e" -fo , 

Jk(0)^0 


(0)=0 


(A.15) 


that in turn gives 


?( 1)=0 


' k ( 0)=0 


DRDk e- ^0 = 1 . 


(A.16) 
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Let us now consider the path integral 


pc{l)=u 

Z{u,u'-,0)-= DWc , (A.17) 

J c(0)—li' 

with 

5'[c,c]= f dr c{t) [dr + iO) c{t) — cc{l) . (A.18) 

Jo 

We can solve the latter in a similar way as above, i.e. by splitting the paths into backgrounds 

[dr + iO)c{T) = 0, c(0) = u' C(t) = u' (A. 19) 

( — + i0)c(T) = 0, c(l) = u C{t) = u , (A.20) 

and quantum fluctuations k{t) and R{t) with the same boundary conditions as above. By setting 

c(r) =e-®®"(u' + «:(r)) (A.21) 

c(t) = + R{t)) , (A.22) 


one gets 


Z{u, u'\0) ■= e 


e uu 


rB(l)=0 


/ k (0)=0 


DRDk e- Jo = e^' 


where we made use of (A.16I. In summary, 


^c(l)=u , , 

/ DcDc e~Jo ‘i'^o(T)[^r+^e)c(T)+Cc(l) _ 
J c{0)—u' 


,e uu 


(A.23) 


(A.24) 


which we thus identify with (A.4|. Thus, by adding a constant term to the action and considering 


a = 1,..., independent pairs of oscillator operators, we Anally get 

(‘C{l)—U 

/ c(0)— 


nc{l)—U / \ / + \ 

J c(0)—u' 


= e 


e u-u'-\-i9r 


(A.25) 


where N := • c = c^“ is the total occupation number operator for the system of N harmonic 

oscillators. 

If, on the other hand, one were to identify 


pc{l) — U / X 

/ DcDc e-Jo dTC-(r)(^r+^e)c^(r)+c■c(l) ^ ^-|g 
J c(0)—u' 


with 


(c^ • c)s := i (c'l’ • c + c • c^) = • c + y 




u'} 


(A.26) 


(A.27) 
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being the symmetrized product, one would then get 


pc{l)—U / \ 

/ DcDc g-/o^‘^'^2“(T)(a^+iejc„(T)+e-c(i) _ 
J c{Q)—u^ 


,e u-u' — iO^ 


and 


pc{l) — U / \ / + \ 

J c{0)—u' 


_ u-u'+i9{s—^-) 


However, setting the eigenvalues of the occupation numbers to be equal 


(c^ • c)s — s = c'^ • c — r 


N 

s = r+- , 


(A.28) 


(A.29) 


(A.30) 


the right hand sides of (A.251 and (A.291 coincide. 
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